Coherent control of atomic excitation using off-resonant strong few-cycle pulses. 
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We study the dynamics of a two-level system driven by an off-resonant few-cycle pulse which has a 
phase jump <f> at t — to, in contrast to many cycle pulses, under non rotating- wave approximation 
(NRWA). We give a closed form analytical solution for the evolution of the probability amplitude 
|(7 a (t)| for the upper level. Using the appropriate pulse parameters like phase-jump <f>, jump time to, 
pulse width r, frequency v and Rabi frequency f2o the population transfer, after the pulse is gone, 
can be optimized and for the pulse considered here, enhancement of 10 6 — 10 s factor was obtained. 
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Modern pulsed lasers produce bursts of light that are 
both ultra-short and ultra-strong, exhibiting durations 
comparable to those of molecular vibrations, and electric 
fields rivaling those near an atomic nucleus [Tj. Attosec- 
ond lasers, emitting pulses with only a few optical cycles 
per pulse [2], hold the promise of controlling the phase 
difference between the carrier wave and its envelope[3]. 

The interaction between strong, broadband electro- 
magnetic fields and atoms, especially laser radiation that 
is tuned far from resonance, is of current interest. Short 
pulses can excite coherence on high-frequency transitions 
that may be used for efficient generation of XUV radia- 
tion [3H5]. Shaped pulses can enhance transient popula- 
tion of the excited state 7' or create optimal coherence 
in TLS [S]. Recently we have found a new analytical 
solution describing the dynamics of a two-level atom un- 
der the action of laser radiation with an arbitrary pulse 
shape and polarization [9J. Furthermore, we have stud- 
ied two mechanisms of atomic excitation: multi-photon 
excitation, and breaking of adiabaticity [4], and we have 
shown [TU] that the latter can be the more efficient. 

Interaction of such ultrashort pulses with a two-level 
atom under rotating-wave approximation does not give 
us the complete picture since the variation of the atomic 
polarization and population within the optical cycle is 
not slow. Thus we should not neglect the contribution 
of the counter-rotating terms in the Hamiltonian while 
studying few cycle pulses interaction with atomic sys- 
tems [TTHTH] . On the other hand if the fields are not too 
strong and the variation of the atomic polarization and 
population within the optical cycle is slow, RWA appears 
to be a good approximation. 

In this Brief report we studied the interaction of few- 
cycle pulses, in contrast to many cycle pulses P^I21j . 
with two-level system. These pulse have a phase jump 
<f> at t = to. Thus they can be characterized by the pa- 
rameters peak Rabi frequency fio, pulse width r, carrier 
frequency v, phase jump <fi and jump moment to along 
with the pulse envelope ( which we have considered gaus- 
sian for the numerical simulation). We present an ana- 
lytical solution for this problem. Using the appropriate 
characterizing parameters, the population transfer, can 



be optimized and for the pulse considered here, enhance- 
ment of 10 6 — 10 s factor was obtained [see Fig. 5(b)]. 

The equation of motion for the probability amplitudes 
for the states \a) and \b) of a two- level atom (TLA) in- 
teracting with a classical field is given as [52] 
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where hw is the energy difference between two levels, 
p is the atomic dipole moment; E(t) = 8{t)cosvt. 
In the rotating-wave approximation (RWA) we replace 
cos(^t)e ±iLJt ->■ e ±iAt /2 where A = w - v [23 j, is detun- 
ing from resonance. Introducing Q(t) = p£{t)/H, Eq.Q 
reduces to 
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which have an integral of motion |C a | 2 - 
define a function fit) = C a (t)/C b (t), Eq.Q yields the 
following Riccati Equation 



0. 
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Then \C a (t\ = |/(i)l/\A + \f(t)\ 2 . Alternatively, we can 
get a second order linear differential equation for C a (t), 
from Eq.Q 



C a {t) 



C o (t) + i?C o (t)=0. (4) 



In this paper we will work without RWA, hence the Ric- 
cati Eq.Q takes the new form as 

/ + in*(t)cos(vt)e- luJt f 2 - in(t)cos{ist)e lLUt = 0. (5) 
The approximate solution for Eq. ^ , in terms of the tip 



FIG. 1: (Color Online) (a) Two-level atomic system, atomic 
transition frequency uj = uj a — W&, detuning A = u — v and 
Rabi frequency Q(t) = p£(t)/h. (b) Few cycle sine (dashed 
line) and cosine (solid line) pulse with gaussian envelope. 



angle 9 is given as [S] 
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where the tip angle 9(t) has been defined as 

e(t) = [ n^cos^y 



(7) 



From Eq.Q we can obtain \C a (t)\ = \f(t)\/y/l + |/(t)| 2 . 
To see how well the approximate solution works, we have 
plotted the probability amplitude |C a (oo)| for a complex 
pulse shape given by Sl(t) — f2o[sech(ai) + sech(ai — 3)] 
[sec Fig 2] . Numerical simulation (dashed) and analytical 
solution (solid) shown in Figs. 2(a,b) are nearly identical. 



A. Pulse with phase jump 

In this section we will investigate the dynamics of a 
two-level atom subjected to few-cycle pulse with a phase 
jump at an arbitrary time t = to- Let us define the Rabi 
frequency f2(t) for our model as 



n(t) 



fi_(t) if t < t , 

n+(t) if t > * , 
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where fi+(t) = e l *T2_(t) and (j> is the phase jump intro- 
duced to the electromagnetic field at t = to- Equivalently 
the tip angle define by Eq.Q takes the form 



9{t) 



0_(t) if t < to, 
8 + {t) if t > t - 
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From the definition of the Rabi frequency Eq.pJ), we 
can easily see that 6 + = e^0_.Thc time evolution of our 
system is divided into two regimes (-co, to) and (fo,oo). 
In both these regimes, the functional form of the solutions 
remains the same. We can write 



/*(*) 



/_(*) if t < to, 
/+(*) ift>to- 
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FIG. 2: (Color Online) Population left in the upper level \a) 
after applying fl(t) = £lo[sech(at) + sech(crf — 3)] pulse as a 
function of the frequency v/u) obtained by numerical solution 
of Eq. (IT]) (dots) and using our approximate analytical result 
Eq.([6j (solid line). In calculations we take fio = 0.04a; and 
a = 0.075a;. In (a) <j> = while in (b) cj> = n, t = 



Eq.(|6]) is the solution for cf> — for the initial condition 
/(-co) = 0. Using the same initial condition we can 
safely write 



f-(t)=i [ dt'$_(<>xp 2 f C-(t 
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As the functional form of /+(t) and /_(t) are the same, 
we can write 



f+(t)=i / di'$ + (t')exp 



C+{t")dt" 
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where = e^$_(t') and C+(t') = C-(*') The con- 

stant c can be obtained by demanding the continuity of 
/^(t) at t = to which gives 



c = i 
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Population transferred to the level \a) during the inter- 
action is given as |C a (oo)| 2 = |/ (oo)| 2 /(l + |/ (oo)| 2 ). 
In order to study the effect of the phase jump let us 
define a relative change in the amplitude 



U(t) - f(t) 



fit) 
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Using Eq.Q, Eq.([l3|) and Eq.Q we get, 
r<t>{t) 



(16) 

The asymptotic value r c / > (oo) can be obtained by t — > oo 
in Eq.(16). We can easily see from the Eq.(16l, that 





2f t U-(t")dt"~ 




2 j£C- (*")<**"" 





r^(oo) attains its maximum value for (f> = n. 
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FIG. 3: (Color Online) Effect of jump time to. (a) Here we 
have plotted the probability amplitude |C a (oo)| against the 
phase jump (f>. Phase jump is introduced at the peak of the 
gaussian envelope, (b) The symmetric influence on the degree 
of excitation with respect to the position of to- The symmetric 
response is lost for shifted gaussian input pulse (c) and (d). 
For numerical calculations we chose Qq — 0.875a;, v = 0.75a; , 
a = 0.331a; and 7 = 1.25a;. 



B. Effect of Pulse parameters: Numerical Analysis 

In this section we will discuss the effect of the pulse 
parameters like phase jump time to, pulse witdth r, de- 
tuning A and peak Rabi frequency f2o on the degree of 
excitation of the upper level \a). For the computational 
purpose we have considered a Gaussian pulse of the form 
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FIG. 4: (Color Online) Effect of a. (a) Probability am- 
plitude I C a (00) I varies in the range from 10 _J ^ 0.7. (b) 
We have plotted |C a (oo)| against normalized pulse width 7T 
for fixed w, v, fio and three combinations of the phase jump 
4> — 0,7r/2,7r. (c) Shows the temporal evolution for the three 
combinations used in (b). For numerical simulation we chose 
Q = 0.875a), v = 0.75a;, 7 = 1.25a; and a = 0.331a;. 
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FIG. 5: (Color Online) (a) Temporal behavior of \C a \ for 
difference combination of <j>. (b) Plot of |C a (oo)| against v/oj. 
For numerical simulation we chose Qo = 0.875a;, to — 0, 7 = 
1.25a; and a = 0.331a). 



fi(t) = ^oe-"" 4 where a = 2^J\d2/t (t is the FWHM of 
the pulse). 

The main result showing the effect of relative posi- 
tion of to, with respect to the peak of the pulse, on the 
atomic excitation is shown in Fig. 3 and Fig. 7(a) where 
we have shown the dynamics of the two-level atom inter- 
acting with few-cycle pulse with a phase jump. In Fig. 
3(a) we have one such scenario of cp = ir/2. Here the 
phase jump is introduced in the field at the peak of the 
gaussian envelope i.e to = and plotted the probability 
amplitude |C (oo)| against the phase jump <fi. Interest- 
ingly the difference in the maximum and the minimum 
value corresponds to A(f> — n. The symmetric nature of 
the atomic excitation is observed in Fig. 3(b) and the 
contour plot Fig. 7(a). With the shifted Gaussian pulse 
Q = n Q e- a2( - t:tt ^ 2 [see Figs. 3(c) and 3(d)] the symmetry 
is lost. Also the effect of the phase jump becomes signifi- 
cant for to within the FWHM of the pulse and gradually 
decreases when to is close to the tail of the pulse. Identi- 
cal response of the system, for 7< ~ 10, is observed for 
three combinations of the phase jump = 0, 7r/2, 7r. 

While investigating the effect of few-cycle pulses on 
atomic systems, the parameter a plays an important role 
for a given value of the carrier frequency v. It determines 
the number of cycles of the field in the pulse. The main 
results showing the effect of a or the pulse width r is 
given in Fig. 4 and contour plot Fig. 7(b). If we look 
at the inset of Fig. 4(a) we see that the probability am- 
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FIG. 6: (Color Online) (a) Temporal behavior of \C a \ for 
difference combination of cf>. (b) Plot of |C a (oo)| against Qo. 
For numerical simulation in (b), we chose a shifted gaussian 
pulse with t s = 1, v = 0.75a;, to = 0,7 = 1.25a; and a = 
0.331a;. fi = 0.875a; for Fig. (a). 
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FIG. 7: (Color online) Contour plot showing the effect of 
pulse parameters like to,T,v and Qq on the population left in 
the excited states \a) in (a), (b), (c) and (d) respectively. The 
influence of the phase jump time to is symmetric with respect 
to to is shown in (a). Parameters used are f2o = 0.875a;, 
v = 0.75a;, 7 = 1.25a;, , to = and a = 0.331a; as required 
appropriately. For (c) we used a = 0.110a; 

plitude |C a (oo)| varies in the range from 10 -5 ^ 0.7. In 
Fig 4(b) we have used three combination of phase jump 
<j> {(f> = 0, 7r/2, 7r) to study the effect of a on the degree of 
excitation. For lower pulse width (2 < 7T < 15)0 = 7r/2 
creates more excitation than <f> — or ir. 

In order to study the effect of detuning A we have plot- 
ted the response of the system in terms of |C a (oo)| for 



the three combination of <\>. Fig 5(a) shows the tem- 
poral behavior while Fig. 5(b) gives the information 
about steady-state population. The probability ampli- 
tude |C (oo)| varies in the range from 4.4 x 10~ 4 ^ 0.4 
for = and 5 x 10~ 5 - 0.9 for <p = tt. When |C (oo)| 
is — 4.4 x 10~ 4 for 4> = we have |C (oo)| — 1 for <f) = n, 
thus we have an enhancement of 10 6 — 10 8 factor in the 
population transfer by introducing a phase jump of n at 
the peak of the envelope function. 

The effect of the peak Rabi frequency fto on the degree 
of excitation of the upper level in shown in Fig 6 and the 
contour plot Fig 7(d). While Fig 6(a) shows the temporal 
behavior of \C a \ on the other hand Fig 6(b) gives the 
information about the population left in the upper level 
after the pulse is gone. We see that for some choice of fin. 
4> = has the maximum effect while for some (f> = 7r/2 is 
dominant. 

In conclusion, we have studied few-cycle pulses, with a 
phase jump <j> at t = to, interacting with a two- level atom. 
This interaction is investigated without the rotating- wave 
approximation and we present an approximate solution 
for the probability amplitude C a (t) of the upper level. 
The approximate solution not only works well with multi- 
cycle pulse [9] but it is also in excellent agreement for 
few-cycles pulses [see Fig 2] . Using the appropriate pulse 
parameters (/), to, a and f^o the population transfer, after 
the pulse is gone, can be optimized and for the pulse 
considered here, enhancement of 10 6 — 10 8 factor was 
obtained [see Fig. 5(b)] 
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